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I. INTRODUCTION 



The problem of the scalar /o(980) and ao(980) mesons has come into a central problem of 
hadron spectroscopy up to charm. The point is that these states posses peculiar properties 
in terms of the naive {qq) model, see e.g. the reviews 0,0. Whereas, all their challenging 
properties can be understood in the framework of the four-quark {q^if ) MIT-bag model 
Q. Along with the q^(f nature of the ao(980) and /o(980) mesons the possibility of they 
being the KK molecules is discussed [^. Furthermore, probably, the /o(980) and ao(980) 
mesons are eye-witnesses of confinement 0. 

During years of time there was established by efforts of theorists ||^,^ (see also references 
quoted in |^) that the research of the decays (p 7/0 •yTTTc and 709 jrin could 
play a crucial role in the elucidation of the nature of the scalar /o(980) and ao(980) mesons. 

At present, the investigation of the (p —>■ 7/0 —>■ '~f'n"K decay by detector CMD-2 is being 
carried out at an upgraded e"'"e~-collider VEPP-2M in Novosibirsk. Another detector, SND, 
aiming to study the decays under consideration has put into the operation at the same 
facility. Moreover, in the immediate future in Frascati there is expected the start of the 
operation of the (p factory DA$NE which, probably, makes possible studying the scalar 
/o(980) and ao(980) mesons in the exhaustive way. 

It seems clearly that the definition of theoretical predictions prior to data securing is a 
natural prerequisite to clear up a mystery of the scalar mesons. 

In all theoretical papers (see [|| and references quoted there) except for [0], where the 
g^g^ nature of the scalar mesons was considered, the approximation of narrow widths of the 
/o(980) and ao(980) mesons was used for their visible widths are 25 — 50 MeV. Moreover, 
experimenters used this approximation to give upper limits to branching ratios of the decays 
700 and 7/0 0. 

But recently [|l^ there was shown that narrow resonance approximation in this instance 
is not valid and predictions of decay branching ratios gained in narrow resonance approxi- 
mation are at least two times overstated. 
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In this connection, we study the scope for production of scalar molecules in the (f) radiative 
decays beyond the narrow resonance width approximation. 

In Section 2 we introduce the formulae needed for our investigation and discuss correc- 
tions for finite widths to propagators of scalar mesons. 

Section 3 is devoted to the model of the — > 700 (/o) transition amplitude. In Section 3. A 
we state (more precise, make clear) the model of the production of the scalar KK molecules 
in the radiative decays. In Section 3.B we calculate in the analytical form the imaginary 
part of the — >■ 700 (/o) transition amplitude that gives near 90% of the branching ratios 
in the tttt and 7771 channels. In Section 3.C the real part of the (f) 7^0 (/o) transition 
amplitude, dominating in the KK channels, is derived partly in the analytic form partly in 
the integral one suitable for the simulation of the experimental data. 

In Section 4 we give the numerical results of our analysis. We found that in the case of 
the molecule nature of the /o(980) and ao(980) resonances BR((f) — > 7/0(00) ^7nr{nr))) ^ 
(1 - 2) X 10-5 , BR{(P ^ 7(/o + ao) ^K^^K-) < 10"^ and BR{(f) ^ 7(/o + oq) ^ 
^K^K^) < 10~^. The mass spectra for the nn , nrj , K^K~ , K^K^ are presented. The 
phase of the scalar resonance production amplitude, causing the interference patterns in the 
reaction e~^e~ 77r+7r~ in the 4> meson mass region, is calculated. 

Section 5 reviews our results and discusses experimental perspectives. 

II. FUNDAMENTAL PHENOMENOLOGY OF ^ RADIATIVE DECAYS INTO 

SCALAR MESONS 

Let us introduce the amplitudes 

M{(l> ^ 7i?; m) = gR{m){e{cl>) ■ 6(7)) , = ao , /o , (1) 

where e(0) and 6(7) are polarization three- vectors of the 0-meson and the 7-quantum in the 
0-meson rest frame, m is an invariant mass of two pseudoscalar mesons a and b produced 
in the R ^ ah decay. 
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In e+e~ collisions the 0-meson is produced transversely polarized relative to the beam 
axis in the center mass system. That is why the amplitudes from Eq. (1) lead to the angular 
distribution 



W{0) = -(l + cos^e) 



(2) 



if one is not interested in a polarization of the 7-quantum in the reaction e'^e~ — > — > jR. 
In Eq. (2) 9 is the angle between the momentum of the 7-quantum and the axis of the 
beams. 

According to the gauge invariance condition, the decay amplitudes are proportional to 
the electromagnetic (electric) field, that is the photon energy (cu) in a low energy region: 



gnifn) uj X const , 



(3) 



if m — > m^, cu — (l/2)m0(l — m^/m^) — > 0. 

The decay width in narrow scalar resonance width approximation 

r ^ 7i?;m = (1 2 • 

3 Att 2mA, ml 



The physically measurable partial widths 



r(0 ^ 7i? ^ 706) = - / 

TT Jma-\- 



mdm 



ma+rrib 



mV{R — > ah] m)r(0 — > ^R] m) 
\DR{m)\^ 



for ah = 7171, 71^7]. 

For ah ^ K+K-, K^K^ 



(4) 



(5) 



r(0^7(ao + /o)^7i^+i^") = 

- / m2r(/o ^ K+K-; m)r(0 ^ 7/0; m) 

TT J2m^+ 



1 _^ 9ao{m)gaoK+K- 



dm , 



r(0^7(ao + /o)^7i^"i^") 
- / m2r(/o ^ X^^O; m)r((/. ^ 7/0; m) 



+ 



9aQ[m)gaoKOKO 



Dfo (m) gfo {m)gf^KOKO ^ao (^) 



dm , (6) 



where 1/ DR{m) is a propagator of a scalar meson. 
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The width of the decay of the scalar meson R into two pseudoscalar meson state ab with 
an invariant mass m is 

TiR^ab;m) = ^-pabim) , 

lOTT m 

PabiTn-) = — m^/m2)(l — m?_/m'^) , m± = rria ± mi, . (7) 
The final particle identity in the vr^vr^ case is taken into account in the determination of 

gfoirOnO- 

In Eq. (6) it is well to bear in mind isotopic symmetry 

dfoK+K- = g/oKORo ! QaoK+K- = —daoKORo ■ (8) 

The use of narrow scalar resonance width approximation in the case under consideration 
is not valid for two reasons iTDI . 



The first and major reason is connected with soft (according strong interaction standards) 
photons. From Eqs. (3), (4), (5) and (6) it follows that the right slope of the resonance is 
suppressed at least by the factor {u/uq}^, where uq = (1 - M|/m^)/2 and Mr is the 
resonance mass. As was shown [l^ it leads to the suppressions of the integral contribution 



from the right slope of the resonance at least by a factor of five for the decays into tttt , vr?] 
channels and at least by a factor of fifty for the decays into the K^K~ , K^K^ channels. 
So, the physically measurable widths Eq. (5) are caused practically fully by "a half " of 
resonance, that is by its left slope in the channels vrvr and ttt]. Only for this reason the results 
(r(0 '-fR; Mr)) obtained in narrow width approximation are overstated two times. 

The second reason is connected with finite width correction in propagators of the scalar 
mesons. Let us take the generally accepted Breit - Wigner formulae. 

If m > 2mK+ , IrriKo , 

1 _ 1 

DR{m) m\ — m'^ — i{r 0(771) + rj^j^{m))m 

TkkM = ^^^(V^l - A7nl,/7n^ + v'l - 4m2,o/m2)l . (9) 

If 27nK+ < 771 < 2771x0 , 



1 



(10) 



When 2mK+ , 2m/fo > m , 

1 1 



DR(m) ^2 _ ^2 I 



m| - m2 + ""l^^^- {^Am\+/w? - 1 + ^Am\^/m? - 1) - zro(m)m 
Where the width of the decay of the scalar R resonance into the ttt] or vrvr channels ro(m) 

is determined with Eq. (7). 

Since the scalar resonances lie under the KK thresholds the position of the peak in the 

cross section or in the mass spectrum does not coincide with tjir as is easy to see using Eqs. 

(9)-(ll). That is why it is necessary to renormalize the mass in the Breit - Wigner formulae, 

Eqs. (9)-(ll), 



2 



m 



2 _ /\#2 9rk+K 



R 



Ml - '-^^{^^rnl.lMl - 1 + ^Aml^/Ml - 1) , (12) 



where M|, is the physical mass squared ( = 980 MeV and Mf^ = 980 MeV ) while m|, is 
the bare mass squared. So, the physical mass is heavier than the bare one. This circumstance 
is especially important in the case of a strong coupling of the scalar mesons with the KK 
channels as in the four-quark and molecule models. Meanwhile it was not taken into account 
both in fitting of data and theoretical works excluding the papers ||l|J^,p!l|] . It should be noted 
that Eqs. (9)-(ll) are applied only in a resonance region. They, for example, have incorrect 
analytic properties at = 0. The expressions, in which this defect is removed, can be 
found in [^^^ . 



Notice that when the scalar resonances lie between the K^K and K^K^ thresholds it 
is necessary to renormalize the mass in the Breit - Wigner formulae in the following manner 



m 



R 



(13) 



IQlT 

The coupling constants in the molecule model PJ8|,p!^: 



^??oi^+x-/4vr = glK^K-f'i^ = 0.6 GeV . (14) 
Notice that in this model MR-mR = 24(10) MeV for Mr = 980(2m^+) MeV. 
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III. MODEL OF KK SCALAR MOLECULE PRODUCTION IN </> RADIATIVE 

DECAYS 



A. KK loop production of extended scalar mesons 

Unfortunately, at present it is not possible to construct the truly relativistic gauge in- 
variant model in the case of a KK scalar molecule production in the radiative decays 
for the nonrelativistic nature of a wave function of a molecule. But it is possible to "rela- 
tivize" a model constructed in the KK molecule rest frame and to place a gauge invariance 
constraint. 

In the rest frame of a scalar KK molecule we consider the production mechanism de- 
scribed by the diagrams in Fig.l, where the cross in the K~^ K~ fo{ao) vertex points to a 
couphng of the K'^K~ state with a extended meson. 

In the scalar molecule rest frame, p = q, the decay amplitude 

T{p,q) = M{p,q)-M{p,0) , 

M{p, q) = Mi(p, q) + M^ip, q) + M^{p, q) , (15) 



where p and q are four- momenta of the meson and the photon respectively. The amplitudes 

{271)^"^^^ ^' D{k)D{k - p + q) 



M,{p, q) = iegnK.K-g,K.K-e^{ct>)e'^{i) j ^.m)^^^^^^^tj^) ' (l^) 
M2{p,q) = 

\ f ^tuw + q)^{2k + 2q-p)^ 

M3{p,q) = 

N f d^k (2k-2p + q)J2k-p)^ 

+ q)D{k-p) 

correspond to the diagrams Fig. la, Fig.lb, Fig.lc respectively. The function describes 
a momentum distribution of a K'^{K~) meson in a molecule, D{k) — k'^ — m'^ + iO, e'*(0) and 
£'^(7) are polarization four- vectors of the (j) meson and the 7-quantum, e^/47r = a ~ l/^^T, 
gRK+K- is determined by Eq. (7) and q^k+k- is determined by the following manner 



A'+A'-) = -?i!iiJi-,ru 



1-^1 . (19) 



ml 



3 1677 

The subtraction at g = is gauge invariant regularization that always is necessary in the 
field theory. 

The "relativization" 

T(p,g) = (e(0) ■ e{-i))gR{m) = gR{m)[{e{<P) ■ 6(7)) - (e(0)g)(e(7)p)/(M)]|p=.- • (20) 

gives the relativistic and gauge invariant amplitude. Emphasize that in Eq. (20) e{(j)) and 
6(7) are polarization three-vectors of the 0-meson and the 7-quantum in the molecule rest 
frame. 

When 0(|/i;|) = 1 Eqs. (15)-(18) reproduce the well definite gauge invariant field theory 
expression [0. 

Our model is the practically same one as in P|. Nevertheless, it should be noted that 
in Eqs. (4.21)-(4.23) of the paper ^ the function 0(|fc|) is used to describe a momentum 
distribution of a K^{K~) meson in a molecule whereas a K^{K~) meson has the three- 
momentum equal ±(/c — q/2). So, Eqs. (4.21)-(4.23) of are not quite correct and our 
Eqs. (15)-(18) make clear the molecule production model ^. 

Let us integrate Eqs. (16)-(18) over in the molecule rest frame. 

f d^k 1 

f <Pk ^ {kq)"^' 
M2{p;q) = egRK+K-9cpK+K-e{(j))e{-f) / j7r^(Pi\k\) " 



(27r)3^'' " I q 



+ 



Ek (m2 + 2Ekm) {2uEk + 2kq) Ek {m? - 2Ekm) (2po^fc - + 2kq) 



^In contrast to |8| we treat with an amplitude rather than a S-matrix element. The difference is 
the factor i. 

Besides, note that in Eq. (4.24) of the sign (-) is lost! Probably, it is a misprint otherwise the 
gauge invariance constraint, see Eq. (3), is destroyed. 
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Ek+q{2uj^ + 2ujEk+g + 2kq) {pi + uo^ + 2poEk+g + 2kq) J 

f d^k ^ (kp)'^\ 
M^{p]q) = -egRK+K-g4>K+K-e{4>)?{l) J (^'^d^D I ~ } ^ 

( 1 1 

X — + — + 

V Ek {m? + 2Ekm) {2poEk + + 2kp} E^ {m? - 2Ekm) {2ujEk + 2kp) 

+ 1 ^\ , (23) 

Eu-,{u^+p'-2u:Eu-,-2kp){pl + jP-2p^Ek-g-2kp)) 



where 



P = ^ , Po-qo = 'm , go = |g] = u; = (m, - m )/2m , Ek = yk"^ + mj^+ - iO 



Ek+g = ^{k + q)^ + ml+-tO , Ek-q = ^ {k - q)^ + m%+ - tO . (24) 

It is easily seen that the integration over the angle between k and q (or p) is performed 
analytically. 

B. Imaginary part of amplitude of — > 7/0(00) decay 
Following HJl^ we use 



where /i = 141 MeV. 

The function 0(|A;|) suppresses the contribution of virtual K^K^ pairs . That is why it 
is natural to expect that the imaginary part of the — 7/0(00) amplitude is essential in 
comparison with the real one. The numerical analyses in Section 4 supports this hope. 

So, the imaginary part of the (j) ~^ if o{(^o) amplitude 

ImT{p, q) = Im{M{p, q) — M{p, 0)) = e{(f))e{'y)Imgji{m) = 

= e{(l))e{-f)Im{gR{m) - gR^m^)) (26) 

merits the individual consideration here. Eq. (25) makes possible calculating the imaginary 
part in the analytical form. This calculation is too cumbersome to present it here. Because 
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of this we shall only explain the genesis of the imaginary part of the amplitude and shall 
adduce its analytical form. 

The contributors of the imaginary part of the — > 7/0(00) amplitude are the regions 
of the integration where the denominators in Eqs. (21)-(23) vanish. We treat these zeros 
using an infinitesimal negative imaginary addition to the mass square (m|-+ — iO), see 
Eq. (24). 

When m < 2mK+ 

Imgnim) = ^egnK^K-g.K.K- _ 4,.). |^ (^1- - a){E, + a) ^ 

EiE2m^{12a'^ - m^) - a^m\{rr? + Aa^) \ 



mu} El - 0? muj rn^ + - 4m|^+ 

mlim' - Aa'){E, - E,) 32(m| - Aml^fl^m'^ - Aa^^f \ 

2o?uj^ ?>m^{ml-Aa?f J ^ ' 

where = m],+ -iJ,^,E^ = f- ^Jml-Am],+ , and = f + ^Jml-Am\+ . 
When m > 2mK+ 

r , ^ /^^ 1 f^'A iE,-a){E2 + a) 

ImgRim) = negRK+K-94>K+K- .^t 2 — 7^1 ~Y 7f — V ^ 

[2tc)^ [m^ — Aa^y yoj-^ \ [E2 — a)[Ei + a) 

EiE2m''{12a'' - m^) - a^m\{m'' + Aa^) \ 

Aa^m^ I 



Aml+^ El-a^ 8ml+ ^ - J m'^ - Amj,+ 

-| — — in — — — In ' -)- 

mw El - a? mw + - Amj^+ 

m2(m2 - Aa^){Ei - E2) 32(m2 - Am],+ f'\m' - Aa'f 
2(^ixr' 3m^(m| — 4a^)^ 

uw? um m + ^Jm^ - Am]^+ J 

Notice that 

8 fi'^iml - Amj^+ 



ImgRim^) = egRK+K-94,K+K- 



3n {ml - 4a2)3 ^ 
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C. Real part of amplitude of cj) ^ 7/0(^0) decay 



As will be seen from Section 4 the real part of — > 7/0(00) amplitude 

ReT{p, q) = Re{M{p, q) - M{p, 0)) = e{<j))e{j)RegR{m) = 
= e{(j))e{-f)Re{gR{m) - gR{m^)) 

dominates in the suppressed K^K~ and K^K^ channels. 
When m < 2mK+ 

RegR{m) = egRK+K-g^K+K-j^^ x 



^ I ? — [— arctan — — — — arctan ^ ] + 



ujm{'w? — Ao?Y a ji m ^^Jim 

^TYIj iy'-\- r / (X»o , TTt \ 9n 

■[(^'^ctan— ) — (arctan ) J + 



K+ 



m 



2 



ujm{w? — 4a^)^ ■^Am 



K+-nn? 



2 + 7^ + \ 9 ^) (1 ^ arctan - ) + — — - + 



(m? — ^0^)0? 2uj 2uj'^m jia ji ujm{m? — 4a^) 

Jo rfc2 + u2122a;3l 



(P + /x2)22u;3 y 

uj'^ml+ ^^{Ek + u; + E_){Ek+u;-E+) 



Ek{2Ekm - w?) {Ek + u; + E+){Ek + uj- E_ 

iu^ml+ {E^ - Ek + Lo){Ek - LJ + E+) 



In 77^ TT^ ^ + 



Ekm{2Ek + m) {E+ - Ek + to) {E^ + Ek - uj) 
ml{E, + E,){Ek + E2) {Ek + P0 + E.f 
Ekm{2Ek + m) ^ {Ek + po + E+f 
, mliEk - E,){Ek E,) (Ek -po + E^f \ 1 „ _ , , 



Ekm{2Ek-m) {E^ - Po + E_y ^ j 

where E+ — ^JhP' + a;^ + 2\k\uj + m\+ , E_ — ^Jk'^ + a;^ _ 2\k\uj + m\^+ and 

RegR(m^) = egRK+K-g4,K+K-j^-j^^[^ arctan ^ - 
^ml- Am\+ + ^m^ - Am\+ U{m% - 3mlml+ - 4m^+a^) 

I n : o H 0^2/^2 /l„2^3 / 



- ^ml-Aml+ Sm^ (m^ - Ao?f 

arctan(°)( ^"^^^ + 8/xm|-+ 16m|+/i _^ 

> >a3(m^-4a2) 3a3m^(m^ - 4a2) 3a3(m^ - 40^)2 



11 



When m > 2mK+ 



RegR[m) = egRK+K-9<i>K+K-j^^ x 



^ - [— arctan In ^ J + 



I [(arctan -)2 + - In^ ^ + V^' " _ + 

u;m(m2-4a2)2L^ 4 m - ^m^ - 4m^+ 4^ 



2^ ^2 



[m^-Aa^)a^ 2uj 2uj^m //a // 

2 / 4 . m , mi^+ — a 



2 



ujm{m? — A:a?) 2ujm\{m? — A:a?Y 2a mK+ -\- a- {2m\+ — niy 

1 mK+-a mmK+ + 2a^\^ m+ Jw? - Am\+ 
[773!^- ^ + ] In 



2(m2 - Ao?) 2a^ mK+ + a o^jj? ) ^ _ ^7^2 _ 4^ 



+ /x2)22a;3 \^ ^ Ek{m^-2Ekm) 

, (£;fe + a; + £;_)(£;fe + 6^-E+) 
X m 



2^2 

+ 

X 



K 



In 



(-E- — El- + lv)(Ei- — lV + E^ 



Ekm{2Ek + m) {E+ - Ek + u;){E^ + Ek - oj) 
ml{Ej, + E^){Ek + E^) ^^ {Ek+p^^E_f 

Ekm{2Ek + m) (E^ + po + 

m^(^fc - E^){Ek - E^) {E, - po + E+f 
Ekim? - 2Ekm) {E^ - Po + E_f 

Ek[m - ZEkTu) -ui + y/m^ - Am\+ J J 

It is easy to verify that the integrand in Eq. (33) is nonsingular at 2Ek = m for the 
principal value of the integral was calculated analytically. That is why Eq. (33) is suitable 
for a simulation of experimental data. 
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IV. QUANTITATIVE RESULTS 



Below we consider two variants: i) ro(AfR) = 50 MeV that corresponds to the visible 
width ~ 25 MeV and the partial intensity of the decay into the KK channels BR{R — >■ 
KK) ^ 0.35 in the molecule model, Eq. (14), and ii) Tq{Mr) = 100 MeV that corresponds 
to the visible width ~ 75 MeV and BR{R — > KK) ^ 0.3 in the molecule model, Eq. (14). 
Please see the definition of ro(m) in Eqs. (9)-(ll), (7). 

We present BR{(j) —>■ '~fR;'m) = r{(f) '~fR;m)/r{(j)), where r((/)) is the (j) meson full 
width, in Fig. 2, the phase of —gR{m) ( 5 = axctaii {I mgji{m) / Reg ^{m)) ) in Fig. 3, the 
spectra 

dBR{(f) ^ -fR ^ -fab) dr(0 7/? -> 706) 2 m^r{R ab;m)r{(f) ^ ^R]m) 



(34) 



dm T((f))dm n T((f))\Dji(m)\'^ 

for a6 = TTTT and Trr] for the different TqIMr) in Fig.4 and Fig. 5 respectively. The spectra 



dBR{(P 7(ao + /o) -fK+K-)/dm = dT{ 



r(0 ifo;m) 

m 



7(ao + /o) ^ iK+K-)/r{<f))dm 
2 



Dfoim) 



+ 



Daoim) 



dm 



(35) 



and 



dBR{(f) 7(ao + /o) -fK''K")/dm = dT{(j) -> 7(00 + /o) ^K''K")/T{(f))dm 



-mT(/o ^ K''K";m 

71 



r(0 lfo]m) 

m 



Df^{m) 



Daoim) 



dm 



(36) 



are presented in Fig. 6 and Fig. 7 respectively. 

As with the four-quark model |0 in the molecule case the interference is constructive in 
the K^K~ channel and destructive in the K^K^ channel, see Eqs. (35) and (36), for the 
K^K~ loop diagram production model. 



2 Notice that dr((/) ^ 7/? ^ 77r7r)/(im = dT{(l) ^ -fR ^ 77r+7r+)/(im + dr{(l) ^ -/R ^ 
^■K^-K^)/dm = 1.5 • dT{(f) — > 7i? 77r+7r+)/dm. 
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We study the dependence of the branching ratios under consideration on the resonance 
mass. Our resuhs are hsted at Tables I-III. Calculating Table I we take into account only 
the imaginary part of the scalar resonance production amplitudes. 



V. CONCLUSION 

At = 980 MeV and ro(m/J = 50(100) MeV BR{^ 7/0 ^ -fmr = 1.9(2.2) ■ 10"^ 
At niao = 980 MeV and ro(m„J = 50(100) MeV BR{<p 700 7777/ = 1.9(2) ■ 10"^ 

As seen from Tables 5i?(0 7/o(ao) ^ 77r7r(7rr/)) ^ (1 ^ 2) x 10"^ , 5i?(0 j{fo + 
ao) ^ jK+K-) < 10-6 and 5i?(0 ^ 7(/o + Oq) ^ 7^°^°) ^ lO"'. 

When the branching ratios from Table I is compared with the ones in the narrow res- 
onance approximation on Fig. 2, it is apparent that these latter are at least two times 
overstated. Notice that when the parameter in the momentum distribution 0(|A;|), see Eq. 
(25), fi = 100 MeV the branching ratios decreases three times (the fi^ low is natural for a 
three-dimensional integral if a integrand has a range of the order of fi). 

The imaginary part of the 7/0(^^0) amplitude gives 90% of the branching ratios in 
the 7r7r(7r?7) channel, see also Figs. 4 and 5. 

The real part of the — > 7/0(00) amplitude dominates in the KK channel branching 
ratios , see also Fig.6 (a contribution of the imaginary part of the amplitude in the K^K^ 
channel is negligible). 

During years of time it was believed by some that the decay 'yKsKs might prove an 
obstacle to research on the violation of CP invariance in the decay K^Ks despite the 
fact that it was shown even in 1987 that BR{(j) — > 7(/o+ao) iKsKs) — 6.5-10"^ in the 
case which is the worst case from the standpoint of research on the violation of CP invariance, 
i.e., the case in which the /o(980) and ao(980) resonances are of a four-quark nature {q^(f)- 
An upper limit in a Pickwickian sense BR{(j) 7(/o + cto) iKsKs) — 3.6- 10^'' was given 
r3[| . The following investigations, see, for example, |T^, only confirm these results. Here, 



m 



see Table III, we also confirm that the branching ratio for the decay 'jKsKs is low. 
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The phase of gR{m) is the very important characteristic of the model production because 
it causes the interference patterns in the reaction e^e^ 77r+7r^ in the (p meson mass 
region. Emphasize that just these interference patterns are used to identify the — > 7/0 
decay. As seen from Fig. 3 the phases in the molecule case (the solid curve) and in the 
four-quark (point-like) case (the dashed curve) differ considerably. The calculation of the 
interference patterns under discussion is a rather complex problem that is a valid one for 
further investigation. 
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FIGURES 

FIG. 1. Diagrams for model. We send you immediatly by request to gubin@math.nsk.su in 
*.bmp format. 

FIG. 2. Branching ratio BR((f> — > ^R;m). The doted Hne is the real part contribution. The 
dashed line is imaginary part contribution. The solid line is the full branching ratio. 

FIG. 3. The phase 5 = avctai,n{ImgR{m)/Regii{m)). The solid line presents the m-dependence 
in the KK molecular case and the dashed one in the q^q^ case. 

FIG. 4. Mass spectra of tttt channel. The solid line is the full spectra. The dashed and 
doted lines are the imaginary and real part contributions respectively, a) Tq = 50 MeV. The full 
branching ratio is 1.9 • 10~^. The branching ratio with account only imaginary part is 1.7 ■ 10~^. 
b) Fq = 100 MeV. The full branching ratio is 2.2 • 10~^. The branching ratio with account only 
imaginary part is 2.0 • 10~^. 

FIG. 5. Mass spectra of nrj channel. The solid line is the full spectra. The dashed and 
doted lines are the imaginary and real part contributions respectively, a) Fq = 50 MeV. The full 
branching ratio is 1.9 ■ 10~^. The branching ratio with account only imaginary part is 1.6 • 10""^. 
b) Fq = 100 MeV. The full branching ratio is 2.0 • 10~^. The branching ratio with account only 
imaginary part is 1.75 • 10~^. 

FIG. 6. Mass spectra of K~^K~ channel. Fq = 50 MeV. The solid line is the full spectra. 
The dashed and doted lines are the imaginary and real part contributions respectively. The full 
branching ratio is 1.0 ■ 10~^. The branching ratio with account only imaginary part is 0.9 • 10~^. 

FIG. 7. Mass spectra of K^K^ channel. F^^ = 100 MeV and F/^ = 50 MeV (see Table III). 
The ontribution of imaginary part is negligible and not presented here. The branching ratio is 
2 • 10-9. 



17 



TABLES 

TABLE I. BR{^ iMfo) ^ 77rr?(7r7r)) 



BR * 10^ 


— > 700 — >■ 77r77 




m„ • GeV 


r„ =50- MeV 


r„ = 100- MeV 


TTI.f ' CtPV 


Ff =50- MeV 


Tf = 100- MeV 


0.960 


2.232 


2.298 


0.960 


2.357 


2.564 


0.965 


2.104 


2.176 


0.965 


2.226 


2.435 


0.970 


1.959 


2.0-15 


0.970 


2.077 


2.297 


0.975 


1.788 


1.904 


0.975 


1.902 


2.149 


0.980 


1.575 


1.753 


0.980 


1.685 


1.991 


0.985 


1.292 


1.590 


0.985 


1.396 


1.822 


0.990 


0.919 


1.049 


0.990 


1.017 


1.224 


0.995 


0.697 


0.902 


0.995 


0.792 


1.072 
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TABLE II. BR{(t)^jK+K-) 



ruf^; GeV 


rriao; GeV 


T/o; MeV 


Tao; MeV 


BR * 10^ 


0.980 


0.980 


50 


50 


1.01 


0.980 


0.980 


100 


100 


0.48 


0.980 


0.980 


50 


100 


0.715 


0.990 


0.990 


100 


100 


0.6 


0.980 


0.990 


50 


100 


0.72 


0.990 


0.980 


50 


100 


1.0 


0.990 


0.990 


50 


100 


1.12 


0.990 


0.990 


50 


50 


1.81 



TABLE III. BR{(i)^^K^K^) 



mf^; GeV 


niao; GeV 


T/o; MeV 


r„o; MeV 


BR * 10^ 


0.980 


0.980 


50 


50 


1.57 • 10-3 


0.980 


0.980 


100 


100 


2.127-19-3 


0.980 


0.980 


50 


100 


2 


0.990 


0.990 


100 


100 


9.48 • 10-^ 


0.980 


0.990 


50 


100 


5.1 


0.990 


0.980 


50 


100 


6.5 


0.990 


0.990 


50 


100 


4.3 


0.990 


0.990 


50 


50 


1.1 • 10-3 
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